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摘 要 作物 根系 生长 不 仅 取 决 于 生理 因素 , 还 取决 于 生态 环境 因素 , 而 土壤 水 分 环境 与 作物 根系 生长 之 间 
的 关系 则 是 局 部 灌溉 技术 设计 的 理论 依据 之 一 。 ee 步 探索 影响 棉花 根系 生长 的 主要 因素 ,本 文 在 根 - 冠 
水 量 平衡 的 基础 上 ,结合 作物 系数 与 叶 面 积 的 关系 模型 、 根 长 密度 分 布 函 数 以 及 根系 吸水 效率 函数 ,应 用 动态 
规划 理论 , 建立 了 棉花 根系 生长 模型 , 并 以 桶 栽 棉 花 试验 结果 进行 验证 。 结果 表 明 , 该 模型 纳入 了 土壤 水 分 环 
境 、 大 气 蒸腾 力 和 叶 面 积 等 影响 根系 生长 的 因子 ， 具 有 揭示 根系 生长 耗 水 机 理 的 作用 。 该 模型 模拟 出 的 棉花 
总 根 长 变化 趋势 与 实测 结果 基本 一 致 ， 当 以 多 年 月 平均 参考 蒸 散 量 (ETo) 作 为 输入 条 件 时 ， 模拟 结果 总 体 误 差 
为 15.41%， 可 以 用 于 工程 设计 。 对 模型 敏感 性 分 析 结 果 表 明 ， 所 建 模型 能 够 反映 棉花 根 - 叶 生长 的 同步 性 ，b 
及 进入 生殖 生长 期 以 后 根 一 叶 之 间 的 水 量 平衡 关系 。 棉 花 根 系 生长 对 土壤 水 分 环境 变化 的 敏感 性 高 于 对 叶 面 
职 变 化 的 敏感 性 ,体现 了 棉花 根系 生长 的 机 理 ， 建 模 方法 可 行 . 本 文 研究 成 果 对 完善 局 部 灌溉 技术 中 灌溉 制度 
的 设计 理论 具有 重要 意义 。 
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Simulation method of cotton root length growth based on dynamic 
programming theory” 
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Abstract Crop root length and soil moisture distribution are important determinants of crop root water uptake potential. Crop 
root length changes with growth stage, which also requires changes in soil moisture environment. Therefore, establishing a 
root growth model to predict crop root growth conditions under normal water use has theoretical and application significance 
regarding the determination of irrigation quota and soil moisture environment indicators, which are required to design 
appropriate local irrigation technology. However, current root growth simulating models are more suitable for the determina- 
tion of irrigation quotas of whole irrigation technology. These models, which are mostly statistical, can not completely reflect 
the relational mechanisms of the growth of plant root system, crop water use and soil moisture environment. In view of the 
above and, a root-canopy water balance combined with crop coefficient vs. leaf area relationship model, root length density 
distribution function and root water uptake efficiency function, a cotton root growth model was developed based on dynamic 
programming theory and verified by experimental results of a barrel-cultivated cotton. The main results of the study showed 
that the model well accounted for the effects of root growth factors such as soil moisture environment, atmospheric 
transpiration rate and leaf area, which revealed to a certain extent the mechanism of crop water use due to root growth. The 
growth characteristics of cotton root length simulated by the model were consistent with measured dates in the barrel 


experiment. When multi-year average monthly mean reference evapotranspiration (ETo) was used as input condition, the 


* 国家 自然 科学 基金 项 目 (51069012) 资 助 

** 通讯 作者 : 李 明 思 ,主要 从 事 灌溉 原理 与 技术 应 用 研究 。E-mail: leemince@sina.com 
陈绍 民 ， 主 要 从 事 灌溉 原理 与 技术 应 用 研究 。E-mail: shaomin_ly24@126.com 
收 稿 日 期 : 2015-07-20 ”接受 日 期 : 2015-11-05 

* This work was supported by the National Natural Science Foundation of China (No. 51069012). 


** Corresponding author, E-mail: leemince@sina.com 
Received Jul. 20, 2015; accepted Nov. 5, 2015 


http:/www.ecoagri.ac.cn 


第 1 期 


陈绍 民 等 : 基于 动态 规划 理论 的 棉花 根 长 生长 模拟 方法 ChinaXiv 合 作 期 利 


overall error of the simulation result was 15.41%. Therefore the model was applicable in engineering designs. Based on 


sensitivity analysis, the established model well reflected the synchronization between the growth of cotton root and leaf, as 


well as the water balance between root and leaf after entering the reproductive period. The sensitivity of cotton root growth to 


changes in soil moisture environment was higher than to changes in leaf area, reflecting the processes of cotton root growth 


and the feasibility of the modeling method. The research significantly improved the design theory of irrigation systems for the 


development of localized irrigation technology. 
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作物 根系 的 分 布 和 生长 过 程 是 作物 生理 生态 领 
域 、 作 物 栽培 和 灌溉 领域 较为 关注 的 问题 ， 尤 其 是 
在 局 部 灌溉 技术 设计 中 ,作物 的 根系 分 布 和 生长 过 
程 是 确定 土壤 有 效 湿 润 区 的 依据 。 所 以 ， 长 期 以 来 
学 者 对 作物 根系 生长 过 程 开 展 了 大 量 研究 工作 (“1。 
众多 成 果 表 明 ， 随 着 生育 阶段 的 演进 ， 大 部 分 作物 
的 根 重 呈 现 “S” 形 生长 特点 ， 即 生育 初期 根 重 增长 较 
慢 ， 然 后 快速 增长 ， 在 中 后 期 达到 最 大 值 ， 后 期 根 
重 减 小 , 所 以 根 重 随 生 育 期 变化 的 整个 过 程 类 似 于 
Gaussian 分 布 ; 也 有 一 些 作 物 生育 后 期 根 重 不 减 小 ， 
而 是 增长 十 分 缓慢 ,整个 生长 过 程 符合 Logistic 模 
型 0 1。 根 长 随 生 育 阶 段 的 生长 过 程 也 基本 上 呈现 
“S” 形 特点 ("小 根 长 生长 速率 为 抛物 线形 后 1]。 也 有 
学 者 对 棉花 (Gossypium spp.) 的 根系 生长 过 程 进行 了 
研究 上 ,结果 显示 棉花 的 根 长 生长 和 根 重 生长 也 都 
呈现 “S” 形 特点 ,而且 在 土壤 水 分 环境 较 好 的 条 件 下 ， 
根 长 生长 过 程 基 本 符合 Logistic 模型 ; 而 存在 土壤 
水 分 胁迫 时 , 根 长 生长 近似 Gaussian 分 布 。 这 些 研究 成 
果 对 于 认识 作物 根系 生长 规律 有 极 大 的 帮助 ， 在 此 成 果 
基础 上 人 们 对 根系 生长 模型 开展 了 大 量 研究 "9, 
较 多 的 是 建立 根系 治 土 层 深度 的 生长 模型 "中, 发 
现 根 长 密度 治 土 层 深 度 方向 呈 负 指数 函数 减 小 ， 即 
地 表 处 往往 根 长 密度 较 大 ， 越 向 土 层 深 处 ， 根 长 密 
度 越 小 。 这 一 类 的 模型 发 展 相 对 较为 成 熟 ， 主 要 以 
作物 冠 层 蒸腾 与 根系 吸水 之 间 的 平衡 原理 为 基础 而 
建立 。 目 前 针对 描述 根系 随 生 育 阶 段 生 长 过 程 的 数 
学 模型 的 研究 尚 少 ， 其 中 以 分 形 理论 为 基础 、 采 用 
微分 L- 系 统 对 作物 根系 几何 形态 建 模 并 模拟 根 长 生 
长 过 程 取得 了 很 好 的 效果 % 小 也 有 学 者 中 利用 人 
工 神经 网 络 理论 建立 了 膜 下 滴灌 棉花 根 长 随 生 育 阶 
段 的 生长 模型 ， 得 到 的 棉花 根 长 生长 符合 Logistic 
模型 ， 并 与 试验 结果 相 吻 合 。 但 是 ， 这 些 模拟 作物 根 
长 生长 的 成 果 基 本 都 是 以 统计 学 方法 建立 的 ,大 多 
不 考虑 作物 生长 中 的 生态 环境 作用 ， 也 没有 反映 作 
物 根系 生长 与 耗 水 以 及 土壤 水 分 环境 之 间 的 机 理 关 
系 ， 难 以 根据 环境 因素 对 作物 根系 生长 过 程 进行 预 
测 。 传 统 的 灌溉 制度 理论 在 制定 灌水 定额 时 只 考虑 
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作物 根系 治 土 层 深度 方向 的 分 布 状况 ， 并 未 考虑 不 
同 生 育 期 的 根 长 状况 , 这 在 用 于 局 部 灌溉 技术 中 灌 
水 定额 设计 时 存在 明显 的 理论 缺陷 ; 而 描述 根系 生 
长 的 统计 学 模型 很 少 考虑 根系 耗 水 因素 。 所 以 ， 以 
作物 根 一 冠 水 量 平衡 原理 为 基础 ， 建 立 作 物 根系 生 
长 过 程 模型 ， 对 于 局 部 灌溉 技术 中 灌水 定额 的 制定 
理论 县 有 完善 和 补充 作用 。 本 文采 用 根系 虚拟 研究 
方法 51， 以 作物 冠 层 蒸腾 与 根系 吸水 之 间 的 平衡 原 
理 和 动态 规划 理论 为 基础 ， 建 立 滴 灌 条 件 下 的 棉花 
根 长 随 生 育 阶段 生长 的 模型 ,为 实现 根系 生长 模拟 
研究 提供 参考 。 
1 模型 建立 
1.1 基本 模型 
作物 根系 生长 不 仅 受 土壤 水 分 环境 控制 ， 还 与 
冠 层 水 分 蒸 散 强 度 有 关 053。 学 术 界 所 提出 的 根系 吸 
水 模型 大 都 是 建立 在 冠 层 蒸腾 与 根系 吸水 之 间 的 平 
衡 理 论 上 "233， 并 且 普 通 认 为 根 的 干 物质 增长 速度 
与 冠 层 生长 速度 呈 阶 段 性 比例 关系 023， 这 就 意味 着 
根系 的 生长 与 冠 层 的 生长 具有 一 定 程度 的 同步 性 。 
因此 ， 本 研究 利用 贝尔 曼 (R. E. Bellman) 动 态 规划 理 
论 52， 以 水 量 平衡 原理 为 基础 ， 通 过 冠 层 蒸腾 过 程 
推 求 根 长 生长 过 程 。 根 据 贝 尔 曼 动 态 规划 理论 ,可 
以 认为 根系 生长 过 程 的 每 一 步 都 应 该 遵从 其 吸水 速 
率 与 冠 层 耗 水 强度 相等 的 原则 ， 而 且 生 长 过 程 中 根 
系 会 主动 在 吸水 较 容 易 的 土壤 湿润 区 内 生长 ， 即 次 
级 根 段 主动 寻找 最 优 的 土壤 水 分 环境 、 温 度 环 境 和 
肥料 环境 。 假 定 土壤 温度 环境 和 肥料 环境 没有 空间 
变化 ,根系 生长 只 受 土 壤 水 分 环境 影响 ， 其 基本 模 
型 为 : 
EI(D=9r(D (1) 
式 中 : ET() 为 作物 冠 层 蒸腾 强度 , mm.d  ; 5.() 为 根 
系 吸 水 速率 , mm.d-; 1 为 出 苗 后 天 数 , d。 
ET() 可 以 实测 ， 也 可 以 根据 参考 作物 潜在 腾 发 
量 换 算 : 
ET(D)=K.(D:Ko(t):ETo(D) (2) 
式 中 : K. 为 作物 系数 ， 取 决 于 作物 叶 面 积 ; Ko 为 土壤 
系数 ， 取 决 于 土壤 水 分 状况 ， 由 于 所 模拟 的 土壤 湿 


http:/www.ecoagri.ac.cn 


58 中 国生 态 农业 学 报 2016 


~ 1/ 第 24 位 


润 区 是 有 利于 根系 吸水 的 湿润 区 ,区 内 的 土壤 水 分 
状况 对 根系 的 胁迫 作用 很 小 , 则 Ko 污 1; ETo 为 参考 
作物 潜在 腾 发 量 ， 也 就 是 大 气 蒸发 能 力 , 由 Penman- 
Monteith 公式 估算 , 在 设计 土壤 灌溉 湿润 区 参数 时 ， 
用 到 的 参考 作物 潜在 腾 发 量 应 该 是 多 年 平均 值 ， 或 
者 是 某 一 气象 频率 下 的 对 应 值 。 
状态 转移 方程 : 
ET(t+1)=K.|[LA(i+1)|-ET, (t+1) (3) 


S(t+1)=a(t+1)-{L(O)+RIAL(+D), rl (4) 
AL= A (5) 


式 中 : S.(1+1) 是 下 一 步 对 应 的 根系 吸水 速率 ; LA 为 
作物 叶 面 积 , cm”; cx(D 为 根系 吸水 效率 , 即 单位 根 长 
的 吸水 通 量 , mmm-; VV 为 根系 所 占据 的 土壤 体积 ， 
m’; R[*] 为 下 一 步骤 的 决策 画 数 ; L(D) 为 根 长 ，m:; 
AL(1+1) 为 下 一 步 根 长 增长 量 , m-d '; 1 为 出 苗 后 天 数 , d。 
1.2 ”作物 系数 
在 作物 主要 生育 阶段 不 受 土壤 水 分 胁迫 的 条 件 
下 ， 作 物 系 数 与 叶 面 积 指数 之 间 大 都 近似 呈 线 性 关 
系 , 而 叶 面 积 指 数 的 变化 过 程 类 似 于 “S” 形 曲线 
即 作物 生育 初期 ， 叶 面积 指数 由 小 到 大 变化 相对 组 
慢 ， 然 后 迅速 增 大 ; 营养 生长 旺盛 时 期 叶 面 积 指数 
达 最 大 值 ， 后 期 又 逐渐 减 小 ， 所 以 叶 面 积 随 生育 期 
变化 过 程 可 用 Gaussian 模型 近似 表达 。 穆 彩 芸 等 已 1 
试验 观测 了 膜 下 滴灌 棉花 的 作物 系数 K. 与 单 株 叶 面 
积 指数 LAI 之 间 的 关系 ,发 现 它们 之 间 近 似 呈 窘 画 
数 关系 。 由 于 叶 面 积 指 数 LAI 与 叶 面 积 LA 呈正 相 
关中， 本 文 建立 K. 值 与 单 株 叶 面 积 LA 的 关系 模型 : 
K.= aLA? (6) 
式 中 :LA 为 作物 叶 面 积 , cm”; a、5。 为 模型 相关 系数 。 


王 允 豆 等 中 对 膜 下 滴灌 棉花 叶 面 积 的 观测 显示 ， 


棉 株 现 蕾 前 期 叶 面 积 生 长 速度 较 慢 ， 从 现 蕾 到 开花 
增长 最 快 ， 开 花 后 到 盛 铃 期 叶 面 积 已 达到 最 高 峰 ， 
以 后 由 于 下 部 叶片 的 衰老 和 脱落 ， 叶 面积 逐渐 降 
低 。 由 此 可 知 单 株 棉花 在 整个 生育 期 内 其 叶 面 积 随 
时 间 的 变化 规律 也 符合 Gaussian 曲 线 变 化 特征 ,所 
以 本 文采 用 的 叶 面 积 与 时 间 的 关系 模型 为 : 
人 0) 

式 中 : a、b、c 为 模型 相关 系数 。 
1.3 ” 根 长 密度 分 布 本 数 

设 定 根系 所 占 土 壤 体积 太 为 一 定 值 , 任意 时 刻 根 
长 密度 在 土壤 体积 VV 中 分 布 是 均匀 的 , 根 长 密度 的 变 
化 仅 与 生长 时 间 相 关 。 众 多 试验 结果 表明 P11， 棉 
花 根 长 密度 随 生长 时 间 变 化 的 动态 过 程 呈 “S” 型 曲 


线 。 在 良好 的 土壤 水 分 环境 、 温 度 环境 和 肥料 环境 
中 ,该 “S” 型 曲线 是 Logistic 模型 rH， 如果 棉花 生长 
后 期 (一 般 指 进入 吐 架 期 ) 受 到 水 分 胁迫 ， 该 “S$” 型 曲 
线 则 是 Gaussian 模型 1。 


a 


Logistic 模型 :qi. (1)= (8) 
l+be® 
-1-5) 
Gaussian 模型 :di (1)=ae 2° (9) 
LD)= (DV (10) 


式 中 : qd1.() 为 根 长 密度 , mm 3; 大 为 根系 所 占据 的 土 
壤 体 积 , m’; a、b、c 为 模型 相关 系数 。 
1.4 ”根系 吸水 效率 函数 

作物 根系 的 吸水 效率 取决 于 根系 吸水 量 与 根系 
的 分 布 、 表 面积 大 小 以 及 根 的 透 过 性 之 间 的 关系 1 
这 些 因 素 均 随 根 系 生育 阶段 而 变化 ， 所 以 特定 条 件 
下 可 以 简化 为 作物 耗 水 量 与 根系 长 度 的 比值 。 苗 期 
根系 较 少 ， 冠 层 耗 水 由 有 限 的 根系 提供 ， 所 以 根系 
吸水 效率 高 ， 随 着 根 长 密度 的 增 大 ， 冠 层 耗 水 由 众 
多 根系 提供 ， 所 以 根系 吸水 效率 下 降 ; 生长 后 期 冠 
层 耗 水 量 减 少 , 根系 吸水 效率 大 幅 下 降 。 因 此， 根系 
吸水 效率 与 生长 时 间 的 关系 豆 数 是 一 个 开口 向 下 的 
抛物 线 方程 。 基 于 此 ， 本文 建 立 以 下 根系 吸水 效率 
与 生长 时 间 的 关系 模型 : 

Q(t)=-a:t +b:t+c (11) 

式 中 : wx(b 为 根系 吸水 效率 , mm.m- ;上 为 出 苗 后 天 数 ， 
d; a、b、c 为 模型 相关 系数 。 
1.5 “模型 求解 约束 条 件 


0<LA<LAn (12) 
0<di.< dm (13) 
0.6x0:<0<0. (14) 
ALCR (15) 


式 中 : LA 为 单 株 作物 叶 面 积 最 大 值 , cm ; dm 为 根 
长 密度 最 大 值 , mm-; 0 为 土壤 含水 率 ; 0 为 田间 持 
水 率 ; R 为 任意 实数 。 

1.6 ”模型 求解 初始 条 件 


LA=6 (16) 
V=p (17) 
du=y (18) 
a= (19) 


式 中 : 5 为 模型 求解 叶 面 积 初始 值 ， cm?; 6 为 模型 求 
解 根系 所 占 体积 ， 本 文采 用 桶 体 体积 ，m*; y 为 模型 
求解 根 长 密度 初始 值 , mm-;“ 上 为 模型 求解 根系 吸水 
效率 初始 值 , mmm '。 
1.7 ”模型 求解 流程 

根据 动态 规划 求解 思路 ， 本 文 所 建 模 型 的 求解 
方法 如 图 1 所 示 , 用 MATLAB 仿真 平台 实现 。 
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输入 初始 条 件 和 系统 基本 参数 


Initial conditions and basic parameters input (LA, ETu d,, V, a, a, b, c, n) 


阶段 数 Stage number 
和 二 1 


输入 决策 变量 
Decision variables input AL(t+1) 


计算 状态 参数 State parameters 
LA(1+1), K.(t+1), A(t+1), L(D, a(t+1) 


状态 转移 方程 State transition equation 
ET(t+1), S(t+1) 


输出 Output 
3AL(t+1) 


1 棉花 根 长 生长 模型 程序 求解 框图 
Fig.1 Flow diagram of cotton root length growth model 
performance 


2 ”模型 验证 
2.1 ”试验 材料 

试验 于 2014 年 5 一 9 月 在 石河 子 大 学 水 利 建筑 
工程 学 院 实验 中 心 进行 (86°%03'27' 玉 ，44°18'25”"N, 海 
拔 451 m)。 供 试 棉花 品种 :新 陆 早 48 号 ;, 试验 用 桶 规 
格 为 (顶部 内 径 x 底 部 内 径 x 高 ): 45 cmx35 cmx53 cm 
( 苗 期 和 鞋 期 )、48 cmx37 cmx55 cm( 花 期 、 铃 期 和 吐 
架 期 ), 试验 用 土 为 砂 壤土 , 风干 碾 碎 过 2 mm 筛 ， 分 
层 装 入 桶 内 , 干 容重 1.39 g.cm- 3， 孔 隐 度 46.14%， 
田间 持 水 率 ( 体 积 合 水 率 )26.06%,， 每 桶 定植 1 株 棉 
花 且 表 面 覆 膜 。 试 验 设 计 苗 期 、 蕾 期 、 花 期 各 4 个 
重复 , 但 根系 采样 时 发 现 没 必 要 将 所 有 重复 都 取 
完 。 苗 期 采样 1 个 , 结果 发 现 苗 期 的 根 长 和 根 长 密 
度数 据 都 太 小 ， 在 模型 中 模拟 意义 不 大 ,没有 继续 
采样 。 所 以 根据 文献 [25] 的 研究 经 验 ， 将 所 剩 处 理 重 
复 保 留 至 蕾 期 采样 。 鞋 期 采样 3 个 重复 ,花期 采样 4 
个 重复 ， 铃 期 和 吐 架 期 各 采样 2 个 重复 。 生 育 阶段 的 
具体 划分 及 取样 时 间 如 下 : 苗 期 ， 从 出 苗 到 现 董 阶段 ， 
一 般 经 历 45 d 左 右 ， 试 验 过 程 中 在 出 苗 后 40 d 进行 苗 
期 根系 采样 ; 蕾 期 ， 棉 株 第 一 果枝 叶腋 出 现 三 角形 


花苞 (大 小 3~5 mm) 即 进入 蕾 期 ， 试 验 过 程 中 在 进入 
鞋 期 后 11 d( 盛 鞋 期 )， 即 出 苗 后 56 d 进行 蕾 期 根系 
采样 ; 花期 , 棉 株 基部 任意 一 果枝 上 第 一 灯 花 开放 
即 进入 花期 ,试验 过 程 中 在 进入 花期 后 10 d( 盛 花 
期 )， 即 出 苗 后 76 d 进行 花期 根系 采样 ; 铃 期 ， 棉 株 
第 1 果枝 出 现 第 1 颗 成 铃 (直径 大 于 2 cm) 即 进入 铃 
期 ， 试 验 过 程 中 在 进入 铃 期 后 10 d， 即 出 苗 后 102 d 
进行 铃 期 根系 采样 ; 吐 架 期 , 第 1 果枝 上 的 铃 开 始 
吐 架 即 进入 吐 架 期 ， 试 验 过 程 中 在 进入 吐 架 期 15 d， 
即 出 苗 后 130 d 进行 吐 架 期 根系 采样 。 

桶 栽 棉花 每 次 每 桶 的 灌水 量 为 苗 期 1 L、 鞋 期 2 工 、 
花 铃 期 3 L、 吐 架 期 2L。 为 避免 水 分 快速 治 桶 壁 下 
渗 ， 在 土壤 表层 埋 入 直径 为 30 cm 的 防水 纸 圈 ， 埋 
深 为 5 cm。 降 水 过 程 中 对 棉花 采取 避 雨 措施 ， 保 证 
水 分 处 理 的 有 效 性 。 

2.2 ”测试 方法 

灌水 前 后 采用 烘 干 法 测定 土壤 水 分 垂直 分 布 。 
为 了 减 小 对 根系 的 破坏 ， 在 距 棉 株 15 cm 处 采用 外 
径 2 cm 的 土 钻 治 深度 方向 每 10 cm 取 1 个 土 样 ， 直 
至 桶 底 , 每 次 取 3 个 桶 ， 取 土 后 采用 过 2 mm 得 的 同 
种 土质 风干 土 回填 。 

出 苗 后 10 d 开始 对 试验 棉花 进行 叶 面 积 测量 ， 
每 次 测量 时 间 间 隔 10 d， 其 中 根系 采样 当天 加 测 。 
利用 毫米 刻度 的 直 尺 对 单 片 叶 采用 十 字 交 叉 法 量 取 
长 和 宽 ， 如 此 测量 所 有 展开 的 叶片 (枯黄 以 及 未 展开 
的 叶片 不 在 计算 范围 )， 最 后 按照 公式 F9: LA= 民 系 7 
(LA: 叶 面 积 ; wx: 系数 ; 素 叶 长 ; 7Y: 叶 宽 ) 计 算出 叶 
面积 ， 本 试验 中 计算 采用 的 系数 为 1。 

各 生育 阶段 采用 双向 切片 法 "采样 获得 根系 分 布 
参数 。 取 根 时 以 棉 株 为 中 心 , 按 10 cmx10 cmx10 cm 
的 样 方 在 水 平和 深度 方向 取样 ， 直 至 桶 边 和 桶 底 。 
拣 出 根 样 后 ， 冲 洗 和 干净， 然后 铺 在 有 刻度 的 纸 上 拍 
照 , 用 R2V 软件 将 根系 图 像 矢 量化 , 用 Office Access 
2003 处 理 矢 量化 结果 并 计算 出 根 长 ， 最 后 求 取 每 株 
棉花 的 平均 根 长 密度 。 

2.3 ”结果 与 分 析 

模型 验证 所 用 冠 层 蒸腾 量 根据 参考 作物 潜在 腾 
发 量 ET 换算 ,利用 Penman-Monteith 公式 计算 石河 
子 市 1981 一 2010 年 多 年 月 平均 腾 发 量 ETu 值 如 表 1 
所 示 。 

由 于 叶 面 积 的 生长 过 程 选取 的 是 Gaussian 模型 ， 
根据 前 人 的 研究 成 果 以 及 美国 建立 的 棉花 生长 模拟 
模型 GOSSYM 显示 5， 棉 花 不 同 生育 阶段 的 根 干 
重 与 冠 层 干 重 的 比值 分 别 近 似 为 不 同 的 常数 ， 即 
根 、 冠 生长 具有 一 定 程 度 的 同步 性 ; 另外 , 试验 中 在 
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表 1 石河 子 多 年 月 平均 ET 值 
Table 1 Annual averages of monthly ETo of the experiment area of Shihezi mm'd 


月 份 Month 5 月 May 6 月 June 


7 月 July 8 月 August 9 月 September 


ETo 5.281 6 6.7183 5 


棉花 进入 吐 架 期 进行 了 水 分 胁迫 处 理 ， 所 以 模拟 中 
根 长 密度 di 选用 Gaussian 模型 ， 根系 所 占 体积 三 为 
桶 体 体积 。 由 于 试验 中 所 取 的 棉花 苗 期 根系 太 小 ， 


6.5127 5.565 3 3.453 3 


其 根 长 密度 很 小 所以， 采用 即将 进入 蕾 期 以 及 蕾 
期 以 后 的 数据 进行 比较 。 文 中 其 他 模型 系数 根据 试 
验 所 测 数据 进行 拟 合 , 结果 如 表 2 所 示 。 


表 2 棉花 生长 模型 系数 


Table 2 Model parameters of cotton growth 


模型 类 别 


模型 系数 Model parameter 决定 系数 (R”) 


Model category 部 


Determination 


b 加 coefficient 


作物 系数 Crop coefficient 0.002 961 
叶 面 积 Leaf area 888.789 
根 长 密度 Root length density (Logistic Model) 3 907.329 
根 长 密度 Root length density (Gaussian Model) 4704.357 
根系 吸水 效率 Root water uptake efficiency 0.000 001 61 


模拟 结果 与 试验 结果 相 比 较 如 图 2 所 示 。 


350 


200 


根 长 Root length (m) 


一 © 一 实测 值 Observed value 
一 会 一 模拟 值 Simulated value 


0 30 60 90 


120 
出 苗 后 天 数 Days after emergence (d) 


150 


2 棉花 根 长 实测 值 与 模拟 值 随 生 长 时 间 的 变化 
Fig.2 Observed and Simulated behaviors of cotton root length 
growth 


试验 数据 显示 ,棉花 即将 进入 蕾 期 至 花期 之 间 
(出 苗 后 40~80 qd) 的 根 长 生长 很 快 , 铃 期 (出 苗 后 100 d 
左右 ) 根 长 生长 达到 最 大 值 ， 进 入 吐 架 期 (出 苗 后 110 d 
以 后 ) 根 长 开始 衰减 ,棉花 从 蕾 期 到 吐 架 期 的 根 长 生 
长 过 程 类 似 抛 物 线 变 化 (如 图 2 所 示 )。 模型 模拟 出 的 
棉花 根 长 生长 趋势 与 实测 值 基本 一 致 ， 只 是 在 出 苗 
后 42~83 d( 蕾 期 和 花期 ) 和 出 苗 120 d 以 后 ( 吐 架 期 ) 
两 个 区 间 内 模拟 值 较 实 测 值 偏 小 。 这 是 由 于 模拟 过 
程 中 的 ET 采用 的 是 多 年 月 平均 值 ， 在 棉花 蕾 期 、 花 
期 及 吐 架 期 阶段 其 值 可 能 比 实际 值 仿 小， 使 得 计算 
出 的 叶 面 蒸发 量 小 ， 导 致 根 长 生长 速度 比 实际 的 
慢 。 不 过 在 工程 设计 中 这 样 选取 ET。 是 合理 的 ?1。 
另外 ,模拟 的 棉花 吐 架 期 根 长 生长 衰减 速度 比 实测 


0.925 9 0.934 4 
103.690 36.850 0.989 0 
1 066.710 0.147 0.893 0 
100.684 40.544 0.9250 
0.000 124 0.0316 0.9916 


结果 快 , 是 因为 除了 模拟 过 程 中 ET 的 影响 以 外 ， 
还 由 于 在 模拟 中 设 定 根 长 密度 在 空间 均匀 分 布 ， 当 
棉花 吐 架 期 叶 面 耗 水 减少 时 , 所 需要 的 吸水 根系 数 
量 也 减少 ， 导 致 整个 空间 的 根 长 同步 减少 ,所 以 模 
拟 出 的 总 根 长 衰减 快 。 而 实测 结果 显示 ， 根 系 并 非 
均匀 分 布 ， 根 长 密度 随 土 层 深 度 先 增 大 后 减 小 ; 水 
平方 向 以 根 轴 处 最 大 , 两边 近似 对 称 分 布 ; 根系 密 
集 区 由 于 存在 吸水 竞争 ， 根 长 有 明显 衰减 ; 而 根系 
稀 松 区 的 根 长 衰减 不 明显 ， 所 以 总 体 上 根 长 的 衰减 
晶 度 比 模拟 值 小 。 模 型 反映 出 了 根系 生长 受 土 壤 水 
分 环境 、 大 气 蒸腾 力 和 叶 面 积 的 影响 ,揭示 了 根系 
生长 耗 水 的 机 理 。 
3 ”模型 敏感 性 讨论 

模型 计算 值 与 试验 测定 值 的 相对 误差 随 模拟 进 
程 呈 先 增 大 后 减 小 再 增 大 的 规律 ， 其 中 最 大 相对 误 
差 (26.19%) 出 现在 出 苗 后 56 d( 蕾 期 ， 叶 面积 和 根 长 
密度 生长 速度 加 快 的 阶段 )， 最 小 相对 误差 (4%) 出 现 
在 出 苗 后 102 d( 铃 期 )， 总 体 误差 为 15.41%， 这 是 因 
为 模拟 过 程 中 的 ET 是 多 年 月 平均 值 , 与 试验 年 份 
的 值 有 差异 。 但 是 这 符合 工程 设计 的 应 用 规则 ， 对 
于 用 作 设 计 土壤 灌溉 湿润 区 的 依据 来 说 ， 已 达到 精 
度 要 求 (10%~20%); 该 模型 是 有 效 的 。 

动态 规划 模型 中 状态 变量 叶 面 积 LA 和 根 长 L 
都 是 关于 生长 时 间 t 的 函数 ， 在 同样 的 时 间 变 化 At 
下 ， 叶 面积 的 相对 变化 量 始终 大 于 根 长 的 相对 变化 
量 。 文 献 [29] 认 为 ， 作 物 生 育 初期 以 根系 生长 为 主 ， 
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此 时 根系 生长 速度 比 叶 面 生 长 速度 快 ; 而 进入 生殖 
生长 期 (对 于 棉花 是 指 进 入 蕾 期 ) 则 以 冠 层 生长 为 主 ， 
此 时 叶 面 生长 速度 比 根系 生长 速度 快 ; 生育 后 期 ， 
两 者 的 生长 速度 都 减 小 。 这 表明 本 文 建立 的 棉花 根 
长 生长 模型 比较 符合 棉花 生殖 生长 期 的 根 - 叶 水 量 
平衡 关系 ， 此 时 根系 的 生长 是 对 根 - 叶 之 间 水 分 平 
衡 关 系 变化 而 做 出 的 响应 。 

本 文 动态 规划 模型 共 涉及 到 作物 系数 K.、 叶 面 
积 LA、 根 长 密度 四 和 根系 吸水 效率 a(?)4 个 影响 模 
型 计算 结果 的 状态 变量 ， 但 真正 具有 自由 度 的 变量 
是 LA、di. 和 cx(D)。 以 棉花 桶 栽 试 验 结果 为 对 照 ， 分 
别 对 这 3 个 状态 变量 变动 -15%, 判断 各 状态 变量 对 
模拟 结果 的 影响 效果 ,由 下 式 求 得 : 

p = -15% — La0% 100%| Xo (20) 
Low 

式 中 :p 为 影响 效果 百分率 , %; 为 具有 自由 度 的 状 
态 变量 ; L_1sw 为 状态 变量 变动 -15% 后 模拟 的 总 根 
长 ; Laow 为 未 变动 状态 变量 模拟 的 总 根 长 。 

影响 效果 百分率 计算 结果 表明 : 叶 面 积 LA 和 
根 长 密度 di 在 原 数 值 的 基础 上 变动 -15%， 对 模拟 
结果 的 影响 基本 一 致 (图 3 所 示 )， 都 是 使 根 长 减 小 。 
但 是 叶 面 积 对 根 长 生长 的 影响 程度 (影响 效果 均值 
为 -13.9%) 小 于 根 长 密度 di, 对 根 长 生长 的 影响 程度 
(影响 效果 均值 为 -14.9%); 而 且 在 整个 生育 阶段 ， 
叶 面 积 对 根 长 生长 的 影响 程度 几乎 是 定 值 ; 而 根 长 
密度 对 根 长 生长 的 影响 程度 是 变化 的 ,生育 初期 影 
响 程 度 小 ,生育 后 期 影响 程度 大 ,出 苗 100 d 后 ( 铃 
期 )， 根 长 密度 对 根 长 生长 几乎 是 直接 影响 。 这 反映 
了 3 个 方面 的 问题 : 首先 , 根 长 生长 与 叶 面 积 生 长 
呈正 相关 , 这 与 前 人 研究 得 出 的 “棉花 根 、 冠 生长 同 
步 性 ”结论 是 一 致 的 中 反映 了 棉花 整个 生育 期 叶 
面 蒸腾 与 根系 吸水 之 间 的 生理 平衡 关系 。 其 次 ， 虽 
然 根 长 密度 分 布 和 根 长 生长 之 间 也 呈正 相关 关系 ， 
但 它们 之 间 的 因果 关系 却 是 根 长 生长 影响 着 根 长 密 
度 分 布 。 特 别 是 在 桶 栽 条 件 下 (或 土壤 湿润 空间 有 限 
的 条 件 下 ), 棉花 生育 初期 的 根 长 生长 有 足够 的 空间 ， 
根 长 密度 增加 缓慢 ; 而 棉花 生育 中 、 后 期 的 根 长 生 
和 空间 逐渐 达到 极限 或 接近 定 值 ， 而 为 了 维持 水 量 
平衡 ， 棉 花 根系 二 级 及 三 级 侧根 仍然 要 在 土壤 湿润 
区 中 生长 ,因此 ， 根 长 生长 就 直接 导致 了 根 长 密度 
的 增 大 。 最 后 ， 叶 面积 和 根 长 密度 分 布 对 根 长 生长 


的 影响 程度 反映 出 根系 生长 对 土壤 水 分 环境 更 敏感 


这 符合 根系 生长 的 一 般 规律 中 
在 水 量 平衡 条 件 下 ,单位 长 度 根系 吸水 速率 决 


定 着 维持 水 量 平衡 所 需 的 总 根 长 。 当 根系 吸水 效率 
xD) 在 原 数值 的 基础 上 变动 -15%， 使 模拟 出 的 根 长 
增 大 (影响 效果 均值 为 17.6%)， 而 且 在 棉花 整个 生育 
期 其 影响 效果 几乎 相同 (图 3 所 示 )。 说 明 根 系 吸 水 效 
率 与 根 长 呈 负 相关 关系 , 即 如 果 根 系 吸 水 效率 低 ， 
则 为 了 维持 水 分 平衡 所 需要 的 根 量 就 多 。 


一 一 叶 面积 Leaf area 
一 9 一 根 长 密度 Root length density 
一 太一 根系 吸水 效率 Root water uptake efficiency 


0 20 40 60 80 100 120 140 
出 苗 后 天 数 Days after emergence (d) 
3 棉花 生长 程度 中 具有 自由 度 的 状态 变量 对 模拟 结 
果 的 影响 效果 百分率 
Fig.3 Effect result percentage of degree state variables with 
freedom on simulation results of cotton root length 


比较 3 个 因素 对 棉花 根 长 生长 影响 效果 的 绝对 
值 ， 发 现 根系 吸水 效率 对 棉花 根 长 生长 的 影响 大 于 
冠 层 耗 水 和 根 长 密度 (或 土壤 水 分 环境 ) 对 根 长 生长 
的 影响 , 因此 ,生产 中 采用 合理 施肥 以 及 提高 地 温 
来 提高 根系 吸水 效率 , 促进 作物 根系 及 植株 生长 中 
另外 ， 对 于 棉花 根系 生长 来 说 ,， 叶 面积 和 土壤 水 分 
环境 是 其 生长 的 外 部 影响 因素 ,而 根系 吸水 效率 是 
其 内 部 影响 因素 ， 所 以 它 的 影响 效果 最 大 。 

前 人 建立 的 棉花 根 长 随 生 育 阶段 生长 的 模型 基 
本 是 以 统计 学 模型 为 主 王 ， 并 不 反映 根 - 叶 或 根 - 土 
水 分 环境 的 生理 生态 关系 ,而 叶 面 耗 水 以 及 土壤 水 
分 环境 恰恰 是 制约 根系 生长 的 边界 条 件 。 本 文 基于 
水 量 平衡 原理 ， 用 动态 规划 方法 建立 的 棉花 根 长 生 
长 模型 能 够 将 这 些 制约 条 件 考虑 进去 ,使 模型 能 
反映 一 定 的 机 理性 质 。 尽 管 模型 中 各 状态 变量 仍然 
采用 了 关于 时 间 的 统计 学 醒 数 ,但 是 模型 结构 上 已 
体现 了 叶 面 耗 水 、 土 壤 水 分 状况 、 根 系 空间 等 因素 
对 棉花 根 长 生长 的 作用 。 关 于 棉花 叶 面 生长 过 程 和 
根 长 密度 变化 过 程 的 研究 已 有 大 量 成 果 可 参考 ， 只 
是 关于 根系 吸水 效率 的 研究 成 果 还 需要 丰富 ， 因 此 ， 
该 建 模 方法 对 灌溉 技术 设计 有 帮助 。 


影响 效果 百分率 
Effect results percentage (%) 


4 结论 
1) 本 文 在 棉花 主要 生育 阶段 没有 水 分 胁迫 以 及 
所 有 根系 都 吸水 的 条 件 下 ,以 水 量 平衡 为 原则 ， 利 
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用 动态 规划 理论 建立 了 棉花 根 长 生长 模型 ， 经 桶 栽 
棉花 试验 验证 , 模拟 的 根 长 生长 过 程 与 实测 情况 基 
本 一 致 。 当 以 多 年 月 平均 ETu 作为 输入 条 件 时 ， 模 拟 
结果 总 体 误差 15.41%， 用 于 局 部 灌溉 技术 设计 中 确 
定 土壤 有 效 湿润 区 是 有 效 的 。 

2) 所 建 模型 能 够 反映 棉花 根 - 叶 生长 的 同步 性 ， 
以 及 进入 生殖 生长 期 eA et 
系 。 对 于 影响 根系 生长 的 外 部 因素 来 说 ,模拟 结 
显示 出 , 棉花 根系 生长 对 土壤 水 分 环境 变 和 
性 高 于 对 叶 面 积 变 化 的 敏感 性 , 说 明 本 文 所 建 模 型 
能 反映 棉花 根系 生长 与 冠 层 耗 水 以 及 土壤 水 分 环境 
的 机 理 关系 ， 建 模 方法 是 可 行 的 。 
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